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Abstract: The masses of supermassive black holes in active galactic nuclei (AGN) can be derived spec¬ 
troscopically via virial mass estimators based on selected broad optical/ultraviolet emission lines. These 
estimates commonly use the line width as a proxy for the gas speed and the monochromatic continuum 
luminosity, A L\, as a proxy for the radius of the broad line region. However, if the size of the broad line 
region scales with the bolometric AGN luminosity rather than A L\, mass estimates based on different 
emission lines will show a systematic discrepancy which is a function of the color of the AGN continuum. 
This has actually been observed in mass estimates based on Ha/H/3 and Civ lines, indicating that AGN 
broad line regions indeed scale with bolometric luminosity. Given that this effect seems to have been 
overlooked as yet, currently used single-epoch mass estimates are likely to be biased. 
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1. Introduction 


Supermassive black holes, with masses ranging from 
millions to billions of solar masses, are components of 
probably all galaxies (e.g., Fletcher 2003; Ferrarese & 
Ford 2005). Accretion of gas onto such black holes gives 
rise to the phenomenon of active galactic nuclei (AGN) 
which in turn are able to influence the evolution of their 
host galaxies via “AGN feedback” (e.g., Fabian 2012). 
Accordingly, knowledge of the masses of supermassive 
black holes is of great astrophysical importance. 

For nearby galaxies, black hole masses M can be de¬ 
rived from the dynamics of stars (e.g., Bender et al. 
2005; Oh et al. 2009), of gas (e.g., Walsh et al. 2013), 
or of interstellar masers (e.g., Herrnstein et al. 2005). 
For distant AGN, black hole masses can be derived from 
spectroscopy of broad optical/ultraviolet emission lines. 
Such mass estimates make use of a virial relation of the 
form 




( 1 ) 


where v is the root-mean-squared line-of-sight gas ve¬ 
locity, R is the (effective) radius of the broad line region, 
G is Newton’s constant, and / is a geometry factor of 
order unity. The velocity v is derived from the widths 
of the emission lines. The radius R can be derived 
via reverberation mapping (Peterson 1993). As radii 
derived from reverberation mapping are available only 
for a rather small number of AGN, the monochromatic 
luminosity of the AGN continuum is commonly used as 
a proxy for the radius. The virial mass estimator then 
takes the form 


M = KV 2 (XL x ) a (2) 

where L\ is the specific continuum luminosity at wave¬ 
length A, k is a constant to be calibrated, and a ~ 0.5 


(Kaspi et al. 2007; Bentz et al. 2009). The monochro¬ 
matic luminosities XL\ are estimated from the contin¬ 
uum close to the emission line from which v is derived. 
Estimators of this type, for various lines, have been 
used successfully for deriving black hole masses for large 
samples of AGN (Greene & Ho 2005; Vestergaard & Pe¬ 
terson 2006; Kim et al. 2010; Park et al. 2013; Jun et 
al. 2015). 

Currently used mass estimators are based on the ex¬ 
pectation that the best proxy for the radius of the broad 
line region for a given line is R oc (A L\) a . A priori how¬ 
ever, this is an ad-hoc assumption that requires justi¬ 
fication, and alternative scalings might be possible. In 
the following, I thus explore the hypothesis that the 
sizes of broad line regions do not scale with XL\ but 
rather with bolometric luminosity. 

2. Color Effects 

I begin with the assumption that the correct mass es¬ 
timator based on an optical/ultraviolet emission line i 
is 

Mi = Ki v 2 L a (3) 

where L is the bolometric luminosity of the opti¬ 
cal/ultraviolet AGN continuum and a ss 0.5 for all 
lines. At least in some cases, we can analyze the same 
target using two lines a , b , resulting in two mass es¬ 
timates M a , Mi,. Provided that both mass estimators 
have been calibrated correctly, the results for M a and 
Mb should be identical (within errors) for any given 
target, meaning 


M a _ Kg /UaA 2 
M b K b \v b J 


( 4 ) 


If we replace the bolometric luminosity L by the 
monochromatic luminosity derived from a narrow band 
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centered on a wavelength A located close to the emis¬ 
sion line to be analyzed, A L\, we find a modified mass 
estimator 

= (5) 

with Li = ( XL\)i for convenience and being another 
constant. When analyzing the same target using two 
different lines a, b, we find 

M' a _ n a /-ry\ 2 g a /La \“ 

M' W 6 UJ ' 

Applying Equation (4), we can reduce Equation (6) to 


K_£a/M“ fM a 
K €b\L b ) * o6 UJ 


( 7 ) 


using here £ ob = (£ a /£;>) for convenience. 

Assuming that (i) both mass estimators M aib have 
been calibrated independently and using sufficiently 
representative AGN samples, and (ii) the ratio M' a /M' h 
can be derived for a sufficiently large and representative 
sample of AGN, the ensemble average of the mass ratio, 
will be unity (within errors) by construction: 



This relation implies £ ab = ((L a /L b ) a ) _1 , at least if 
the values M' a IM' h are derived for a sufficiently large 
and representative sample. In case of small samples (as 
frequently the case), £ ob needs to be estimated sepa¬ 
rately from the typical optical/ultraviolet color of the 
AGN continuum. 


M (xv 2 . Likewise, the power law indices of L 1350 and 
L 5100 are consistent with 0.5 within errors. 

From their analysis, Assef et al. (2011) find a sys¬ 
tematic discrepancy between the masses derived from 
the Balmer lines on the one hand and from C IV on the 
other hand. This discrepancy follows the relation (their 
Equation (8)) 


= 10 -v(^) X ( 10 ) 

M CIV V-O5100/ 

where the subscript H denotes either Ha or H/3. De¬ 
pending on which Balmer line and which velocity indi¬ 
cator (line dispersion or line FWHM) is used, best-fit 
values range from 0.51 ± 0.14 to 0.95 ± 0.22 for x and 
-0.11 ± 0.06 to -0.27 ± 0.08 for y. 

If indeed the broad line region radius scales with bolo- 
metric luminosity rather than with XL\, i.e., if Equa¬ 
tion (3) gives the true underlying relation, then Equa¬ 
tions (7) and (10) are equivalent. This can be checked 
in a straightforward manner by identifying either Ha 
or H/3 with “line a” and CIV with “line 6”. Evidently, 
Equation (7) predicts a = x. Indeed, the values ob¬ 
served by Assef et al. (2011) for x are in agreement 
with a ss 0.5 within errors. 

The ensemble averaged optical/ultraviolet continuum 
of quasars is known to follow the relation A L\ oc 
A - O' 56 in the wavelength range 1300-5500 A (Vanden 
Berk et al. 2001). Using a = 0.5, this implies 1 
£, ab « [(5100 A/1350 A )- 0 - 56 ]" 0 5 « 1.45. Compari¬ 
son of Equations (7) and (10) leads to the prediction 
log Cab ~ 0.16 = —y. Within errors, this is in agree¬ 
ment with the values observed by Assef et al. (2011). 

4. Discussion and Conclusions 


3. Comparison to Observations 

Arguably the most careful cross-analysis of different 
mass estimators is provided by Assef et al. (2011) who 
compare mass estimates based on the C IV A1549 line 
to those based on the hydrogen Balmer lines Ha A6565 
and H/3 A4863 (with all wavelengths in A). Assef et al. 
(2011) apply the mass estimators 



drawn from Greene & Ho (2005), Vestergaard & Peter¬ 
son (2006), and Bentz et al. (2009), to a sample of 12 
lensed high-redshift quasars. Here L 1350 and L 5100 de¬ 
note the continuum luminosities A L\ derived at 1350 A 
and 5100 A, respectively; the values of the constants / 
and e depend on if v is derived from the dispersion or 
the FWHM of a line. The relation between M[ ia and 
r>Ha is actually M[ la oc Uy° 6±0 ' 06 , i.e., consistent with 


From the analysis provided in Section 3, it is straight¬ 
forward to see that the observed systematic discrep¬ 
ancy between mass estimates based on different emis¬ 
sion lines is in agreement with AGN broad line regions 
scaling with bolometric rather than monochromatic lu¬ 
minosity. On the one hand, the fact that AGN broad 
line regions are shaped by both the optical/ultraviolet 
and the ionizing Lyman continuum was already noted 
by Peterson (1993) and Baldwin et al. (1995). On the 
other hand, little effort has been made so far to distin¬ 
guish scalings with L from those with A L\ observation- 
ally. It now seems that the data of Assef et al. (2011) 
provide the as yet most direct evidence for a scaling 
of the radii of AGN broad line regions with bolometric 
luminosity. 

Comparisons of black hole masses derived from opti¬ 
cal and ultraviolet lines usually focus on the line widths 
rather than the broad line region radii (cf., e.g., Ho 
et al. 2012 vs. Runnoe et al. 2013); this might ex¬ 
plain why color effects have largely been overlooked so 
far. Even though, an independent hint might have been 
provided by the Baldwin effect (Baldwin 1977), a char¬ 
acteristic anticorrelation between the equivalent widths 

1 I use here {(L a /Lb ) a )~ 1 Rj (L a /Lb )~“. This is possible be¬ 
cause the ratio L a /Lj, is of order unity always and 0.56 a < 1. 
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of ultraviolet emission lines and continuum luminosity: 
whereas equivalent widths are usually discussed as func¬ 
tion of XL\, the actual underlying correlation might be 
the one between equivalent width and bolometric lumi¬ 
nosity (cf., Xu et al. 2008, and references therein). 

When expressed as functions of bolometric AGN lu¬ 
minosity (to ease comparisons of different lines), char¬ 
acteristic broad line region radii can be assigned to in¬ 
dividual emission lines; for the case of H/3 vs. C IV, the 
radii of the former are about three times larger than the 
radii of the latter in general (cf. Chapter 7.1.8 of Net- 
zer 2013). Such a gas stratification is the natural con¬ 
sequence of the differences in ionization potential; lines 
requiring higher excitation or ionization energies have 
to be located closer to the central source of radiation. 
Indeed, studies of the radius- luminosity relationships 
of broad line regions assume (at least implicitly) that 
XL\ is a sufficient proxy for the bolometric continuum 
luminosity (e.g., Kaspi et al. 2007; Bentz et al. 2009). 
Such studies derive radius luminosity relationships, as 
well as black hole mass-luminosity relationships, by 
analyzing radius or mass as function of XL\ and ap¬ 
plying global scaling factors that implicitly include an 
ensemble averaged bolometric correction. However, if 
broad line regions scale with bolometric AGN luminos¬ 
ity, bolometric corrections need to be applied source by 
source. As demonstrated by Assef et al. (2011), ne¬ 
glecting colors and bolometric corrections leads to dif¬ 
ferences up to a factor of four between black hole masses 
estimated from the Balmer lines and C IV, respectively. 
This also suggests that all emission-line based mass es¬ 
timators are biased to varying degrees. 

Assuming a scaling of mass estimates with bolometric 
luminosity, the cross-comparison of non-simultaneous 
virial mass estimates based on different lines is affected 
by the variability of the AGN luminosity. Quasars are 
known to vary in luminosity by factors up to about 
three on time scales of years (e.g., Schramm et al. 1993); 
AGN variability in general obeys red noise statistics, 
meaning that stronger variations occur on longer time 
scales (Park & Trippe 2012, 2014; Kim & Trippe 2013). 
However, such variability cannot introduce systematic 
trends into samples of sources; rather, it increases the 
measurement error for individual sources. In general, 
variability effects are relatively moderate because only 
the square root of the luminosity enters the mass esti¬ 
mate. The individual mass estimates used by Assef et 
al. (2011) come with errors of around 0.3 dex; accord¬ 
ingly, variations in luminosity up to factors of about 
5 would be within the measurement errors (because 
log(V5) = 0.35). 

The discussion provided in this paper, along with 
the one by Assef et al. (2011), illustrates once more 
the fact that any systematic scatter in a given relation, 
like the one between mass estimates based on Ha/H/3 
and C IV, is indicative of hidden parameters. Only a 
careful analysis of the impact of AGN continuum color 
eventually leads to the insight that AGN broad line re¬ 
gions scale with bolometric luminosity rather than XL\. 
(Another recent example is provided by the relation be¬ 


tween Bondi accretion rate and kinetic jet power in ra¬ 
dio galaxies; Trippe 2014.) In turn, blind application 
of “standard” relations despite the presence of obvious 
systematic scatter comes with the risk of producing re¬ 
sults that are seriously biased. 
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